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LOW-THRUST SOLAR ELECTRIC PROPULSION 
NAVIGATION SIMULATION PROGRAM* 

.by ’ 

H. Hagar, Jr. , and ,T. J. Eller 

General . 

Program LOGO is an interplanetary low-thrust, solar electric propulsion 
mission simulation program suitable for basic navigation feasibility studies. 

It employs simple, two-body dynamics to simulate the heliocentric phase (no n- 
body perturbations). Provisions are made to simulate uncertainties in the 
thrust program in a realistic manner, and to assess the effects of these uncer- 
tainties on various navigation strategies . One key feature is the ability to 
configure the dynamic model equations in a number of different ways to account 
for the thrusting uncertainties. Several navigation data types may be simulated: 
Earth-based radar range and range-rate, and on-board celestial observations in- 
volving the sun. Earth, and a specified navigation star. Although gravitational 
perturbations of the Earth acting on the spacecraft are not considered, rotational 
dynamics of Uie Earth. are modeled to account for the significant effects of track- 
ing station motion. 

Several types of information output. are available. These include detailed 
numerical information print output, and both printer plot and digital plot 
features. Further., a limited on-line display capability is available for tele- 
type and CRT terminal facilities. 


■ *This research was supported by the Jet Propulsion Laboratory under Contract 
No. 953147. 



MATHEMATICAL MODELS 


Trajectory Simulation 

The motion of the solar electric propulsion (SEP) spacecraft is assumed 
to ho governed by the gravitational -attraction of the sun (which is assumed to 
bo per foctly known), and the thrust acceleration of the solar electric propul- 
sion system . further, random errors in the. thrust program are assumed to in- 
fluence the spacecraft motion. If only the central force attraction of the sun 
is included, the equations of motion for the SEP . spacecraft are 




+ T 


( 1 ) 


where, as shown in Figure 1, r is a 3-vector of heliocentric position components, 

X, Y,' Z; v is a 3-vector of heliocentric velocity components X, Y, Z; |r| is 

the magnitude of. r; and y is the gravitational parameter of the. sun. T is the 

heliocentric thrust acceleration vector composed of the design thrust acceleration, 

?*, as well as thrust acceleration errors from a number of sources (beam voltage 

and current, grid warpage , deadband control errors, etc.). The heliocentric 

components of T, [T„ I T I T ] , may be expressed in a vehicle centred frame as 

[T IT It], where the two vectors are related by 
x.y.z 5 J 



where ip is. the heliocentric- orientation angle (see Figure 1). 

Random, errors are assumed to occur, in both the magnitude of the- nominal 
thrust acceleration, a*, as well as in the thrust vector pointing angles, y 
and 0 (see Figure ; 1) . Thus* the true thrust acceleration vector, T, differs 
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from the nominal thrust acceleration vector, T* , in both magnitude and orienta- 
tion. The true thrust acceleration, expressed in the vehicle centered reference 
frame (x, y, s) is (see Figure 1) 



where a = a* + 6a. The acceleration error magnitude is simulated by 6a - 

5a sin art + u where 6a and m are constants and where the random variable, 
o a o . 

u , has the statistics E{u } = 0, E{u 2 } = a 2 . For the design mission, T* is 
a a a a 

assumed to be of constant magnitude along the vehicle centered y-axis and the 
nominal values of the pointing angles, y and 6, are assumed to be related as 
shown in Figure 2. a. The radius of the circle is the maximum deviation, sin y * y, 
of the thrust vector from the nominal position; it represents deadband pointing 
errors,. The quantity d = s(t-t^) is the distance of the tip of the thrust 
vector from the point where it last touched the boundary, y. The rate, s, is 
simulated as a constant plus an additive noise component; t is the current mis- 
sion time, and t^ is the time at which the boundary was last encountered. The 
angle, cf), is sampled from a uniform distribution U(-.0708, .866). In addition, 
purely random components are added to s, y, and 0. With these assumptions, 

Eqs. (1) can be integrated to obtain the simulated true trajectory. (Figure 3 
shows the trace of the tip of the thrust acceleration vector as projected onto 
the local x-z plane.) 

Models for Error Compensation 

In the subsequent discussion, it is assumed that the thrust acceleration 
can be separated into modeled and error components, i.e., T = T* + m(t) where 
rn(t) is a 3- vector of thrust acceleration error components, the errors, m(t). 





k+l 



Figure 2.b. Pointing Angles Simulation 
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are approximated by e(t) where e(t) satisfies one of several possible first- 
order or second-order differential equations. The values of e(t) and any un- 
specified parameters in the differential equations which describe e(t) are 
estimated simultaneously with the position and velocity components . Six models 
are available for use as approximations for m(t). See Tables la and lb. 

Model 0 . In this model an arbitrary state noise covariance matrix, Q, 
is added to the differential equation governing the state error covariance. 

The Q-matrix compensates for the process noise and is used to maintain a positive 

definite error covariance matrix. The differential equations for the state vector 

V T r T T-f 
X •= [r v J , are 

r + T* + u (4) 

where u is a random 3-vector with the a priori statistics E{u} = 0, and 

■ t 

E[u(t)u ( t ) ] = q(t)6(t-T).’ This model corresponds to a value of M0D0P = 0 or 
M0D0P = 1. (M0D0P = 1 not fully implemented, do not use ) . . See Table lb for the 
form of state noise covariance logic. 

Model 1 . The thrust acceleration error components are approximated by a 

first-order C-auss-Markov process. The differential equations for- the state 

T - T T T T T T . 
vector, X = l.r v e n ot 3 ], are 

r = v , v — — r + S(e) 

| r | 3 

e + u^ , n = , a = u a , 3 = u g (5) 

where u £ , u^, u^, and u^ are random processes with the following a priori 
statistics 

Elu^} = 0 E{u^(t )u^ ( t ) } = q^(t )6(t-x ) , • e , n 9 a , 3) 





9 - 


(In the program, selection of this model corresponds to a value of the parameter, 
MODOP = 11.) The function, S(e), is defined in Table la. 

Model 2. * This model also employs a first-order Gauss-Markov process to 
approximate the thrust acceleration error components. The form is slightly dif- 
ferent, employing fewer variables (and hence, has less flexibility), to yield the 

T T T T T 

differential equations for the state vector, X = [r v e a ]: 

r = v - , v - - — — r + . S ( e ) 

( 6 ) 

e = -[^a J e + u , a = u 
e 5 a 

where u and u are random with 
e a 

E{u^ } = 0 , E{u^(t )u^(t) } = q^( t )<$ (t-x ) , £e{e,-a} 

(This model corresponds to a value of MODOP = 21.) Again, S(e) is defined in 
Table. la. The matrix is diagonal and its elements are those of the vector, 

a. This model is basically the same as Model 3, except that the appropriate state 
error covariance terms must be set to zero. 

Model 3 . The thrust acceleration error components are approximated by two 
first-order Gauss-Markov processes. The differential equations for the state 
vector, \ =[r v . e n a 3 J, are 

r = v , v = - — — r t S(e,n) 

e = -E"aJ £ + u , d = u 
e a 

(7) . 

n = - n + u , 3 '= u fi 

n » 

where, u , u , u , u 0 are- purely random with 
e n a p 
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E{u £ ) = 0 , E{u £ (t )u £ (t) } = q £ 6(t-T ) , 2,e{£,n,a,g} 

(This model corresponds to MODOP - 21.) The vector function, S(e,n), is defined 
in Table la. - The matrices, F a-J and.F^g-J, are diagonal matrices whose diagonal 
elements form the components of the vectors, a and g. 

Model 4 . The thrust acceleration error components are approximated by a' 

second-order Gauss-Markov process. The differential equations for the state 

T T T T T T T-. 

vector, X=[rvenagJ, are 

r = v , v = - — — r + S(n) 

|r | 3 

e = e - ["3J ri + u £ (8) 


n = e , d - u a , g = u g 

where , and u^ are purely random processes which satisfy the a priori 

statistics 

E{u £ } = 0 , E{u £ (t)uJ(T)} = q £ (t)6(t-x) , £e{e,a,g} 

(This model corresponds to MODOP = 12.)- The vector, S(n.), is defined in 
Table la. F" aj and 6-J are diagonal matrices whose diagonal elements form 
the components of a and g. 

Model 5 . The thrust acceleration error is approximated by a slightly 

different form for the two first-order Gauss-Markov process. The differential 

T T T T T T T 

equations for the state vector X=[rvenotg] are 


r =. v 


v = 


r + S(e,n) 


e = - 


0 0 • 
0 0 


z + u 


a = u 


(9) 


Bj 0 0 

o b| o 

0 0 31 


n + u ? 3 = u„ 

n 3 
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where u , u , u , u D are purely random with 
e r) ql p 

E{u^} = 0 , E{u^(t )u^( x ) } = q^(t)6(t-x) , £e{ e , n ,a , 6} 

(This model corresponds to MODOP = 22.) S(e) is defined in Table la. The 
elements of cu and 8^ are the elements of the vectors, a and 8. 

The Vector Function, S 

In all models, the function, S, may be defined in one of three ways, de- 
pending upon the value of the program parameter, KSW. (See also Figure 1.) 

Table la is a cross-reference matrix showing the form of S for various values 
of MODOP and KSW. (KSW is also called KAX in the program.) 

Observations 

Observations for the orbit determination process consist of range, p; 
range-rage, p; and three celestial angles, A^, defined below. 

A^ = sun-vehicle-planet angle 

A^ = star-vehicle-planet angle 

Ag = sun-vehicle-star angle 

The true and nominal state vectors are used to compute the true and nominal ob- 
servation values. Further, the true observations are corrupted by adding white 
noise to the deterministically computed values. The random components are ob- 
tained by snmpl.itig from known Gaussian distributions. If each observation, 
specified generically as ft, is made at discrete times, t^, the observation-state 
relationship can be expressed as 

a. = G 0 (x . , t. ) + v. do) 

i fii i i 

The observation error, is assumed to have the a priori statistics E{v^_) = 0 , 

E{v n v 0 } = R_ 6. where <5. is the Kronecker delta, 
fii fij ft*. ij ij 
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E arth Ephemeris 

The Earth’s orbit is assumed to be circular, an assumption not unreasonable 
in light of other model dynamics employed. The position and velocity of the 
Earth is determined from the following equations: 

ip = mod (to t , 2tt ) 

Lj 

Xg = C cos ip 
Y r = C sin <p 

( 11 ) 

V = C “e 

X^ = -V sin \p 
Y e 3 V cos ip 

where C is the mean Earth-Sun distance (see Figure 4.). 

Tracking Statio n Motion 

Figure 4 shows the coordinate frames and variables used in defining 
motion of (up to three) tracking stations, due to the Earth’s rotation. The 
Earth is assumed to rotate at a constant rate; no precession or nutation is 
accounted for. From the figure, the tracking station position and velocity in 
heliocentric coordinates are seen to be: 
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where to is the. Earth spin rate, SI is the orbital angular rate of Earth, 

t h 

4k = cot t , and where the subscript, i, represents the i station. 


and 
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-•or tiie geographic rectangular station coordinates x , y , z , D is 
' S i s i s i s i 

' X s + - v s for the ith station i is the right ascension of the i th station 
& i i 1 

at the initial time. The obliquity of the ecliptic is e; in the program 
sin e = .3979 ... and cos e = .9174 ... are coded directly 



Figure 4. Earth Ephemeris and Tracking Station Geometry 
Filter Equations 

The differential equations described by any of the above models may be 
written in the following vector form to obtain the equations of state: 

X = F(X ,t ) , X(t o ) = X o (13) 

The extended form of the Kalman-Bucy filter is used to obtain the state estimate. 

That is, given a previous estimate, an< ^ assoc i a Ted state error covariance 

matrix, P , the estimate and error covariance at time, t 1 , are obtained from the 
k-1 ^ 
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following equations: 


ft. 


X k ' X k-1 + 


R = P, + 

k k-1 


F(X( t ) , x )dx 


'k-1 


P( t )dx 


“k-1 


P = A(t)P + PA T (t) + Q(t ) 


“ T ,T r T T s „T 


K. = P. H [H, P. H + R n , ] 
k k k k k k ftk 


-1 


x k = \ * V a k - W V ] 


V CI -W f k 

where A(t) = 9F(x,t)/3x and = 9G^(x k , t k )/3x, and Q'e(q, q^, q , q 
appropriate non-zero submatrices of the state noise covariance, Q. 


( 14 ) 


a’ V’ the 



Program Description 


M ain Fro gram - LOGO 

The main program provides the controlling logic for driving the overall 
program. All problem input and batch mode control parameters are input here. 
Initialization of program options, and associated logic flags and parameter 
values are set. Plot option logic is set and the plot subroutine is called (on 
problem termination) . Program results are saved on a separately specified output 
file for later use as inputs for continuation of problem at hand. 

Language 


FORTRAN IV 

List of Subroutines used 


FORTRAN IV 


< 


f OUTPUT ^ 
PATH 
UPDATE ' 
MOTION 
ACCEL 
AMATRIX 
HMATRIX . 
OBSERV 
^ GAUSS ) 
MINIPLT 


These are explained in this report 


plotting routine for UT plotting system 


COMPASS 


f RANF \ 
MSAB 
/ MSABT 
MSATB 
^IPRINT^ 


UT2 system library; generates uniform random sequence 
matrix multiplication C = AB 

T 

matrix multiplication C = AB 

T 

matrix multiplication . C = A B 

available for error messages from the above 3 routines 
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Control Cards: 

READPF <number> NAV RMS 

REWIND NAV RMS 

RUN TX I=NAV 

COPYBF RMS LGO 

SETCORE ZERO 

LGO 

Computing System 

CDC 6600, UT-2 Operating System 

Inputs 

Input parameter values are entered under both formatted and NAMELIST input 
(READ) statements. They are explained in Table 2. 

Outputs 

All input parameters are printed immediately after being read. Upon comple- 
tion of problem computations, selected variables are plotted (position, velocity, 
etc.) depending upon the plot option selected. Also, all input parameter values 
are saved on output file TAPE1 , and can be subsequently used, if desired, as the 
input source for a later continuation run. This sequence occurs at the end of 
each run; in order to use TAPE1 as a new input file, it must be saved as a 
permanent file , or punched on cards . 

Unit s 

Any set of engineering units may be employed. All computation and output 
are done in the selected set with one exception. The choice of Earth radii and 
days as units of distance and time results in computation and print output in 
these units; however, plotting is in kilometers and seconds. 

Input File Setup 

Basic problem input data is entered via the NAMELIST INPUTS whose elements 
are listed in Table 2. Since the data is via NAMELIST, any unspecified parameters 
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are automatically assigned values of zero if the SETCORE zero control card has 
been used. Figure 5a shows a sample data deck setup for the case where the a priori 
state error covariance matrix is diagonal. Figure 5b shows the setup for non-zero 
elements in the a priori error covariance. Note, in this case, the addition of 
NAMELIST APCOV, with elements made up of the array, P. Further, the array, STERCOV 
in NAMELIST INPUTS, must be zero in order to read in APCOV. (See Table 3.) 



Figure 5a. Normal Deck Setup 
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On-Line Feature 

The on-line feature of the program is primarily that of presenting a 
running display on an assigned CRT terminal. Selection of the on-line feature 
is done by assigning a negative value to the parameter, PB, in NAMELIST INPUTS. 
Each time the output routine is entered, a colon, :, is displayed indicating 
readiness to receive a keyboard input. The keyboard input is one of the digits, 
0 through 7; their meaning is shown in Table 3. 


K/B Input 

Definition 

0 

Return and continue problem computations 

1 

Print state deviation 

2 

Print state vector 

3 

Print thrust acceleration (local frame) 

4 

Print error covariance trace and RSS 

5 

Print ’ observation error and H-matrix 

6 

Discontinue CRT display after next 0 entered 

7 

Revise problem termination time to the current 
time plus one integration step 


Table 3. On-Line Feature . Keyboard Inputs 

Print Output 

Printed output data consists of true (simulated) and nominal (estimated) 
state vector and state vector deviation, Earth state vector, error covariance 
state noise covariance, filter gain, observations, observation error, and the ob- 
servation matrix, d (observation)/9 (state). A sample block of print is shown in 
Figure 6. 
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Flow Diagram 

Figure 7 is a functional flow diagram of the main program logic. 

Variables 

Tables 2a and 2b define the NAMELIST INPUTS and APCOV parameters. Tables 4 
through 14 define the labeled common blocks (there is no blank COMMON). .Table 15 
defines those internal variables not contained in Tables 2 through 14. These 
internal variables are in alphabetical order. 

Subroutine OUTPUT 

This subroutine provides the logic for controlling print output as well as 
the on-line CRT display logic. Data to be output is communicated to the routine 
via labeled COMMON. Figure 8 is a functional flow diagram of the subroutine. 
Subroutine PATH 

Subroutine PATH performs the numerical integration of the equations of 
motion, and stores the data points for plotting. Two orders of integrators are 
available: Fourth-order Runge-Kutta-Gill, and second-order Runge-Kutta. The 

type of integrator is specified by the "hard-coded" value of the parameter, IGO. 

IGO = 2, Second-Order Runge-Kutta 

Pui = y i + 

y i + i = y i , [f<y i*i- t iti ) - f(y i’V ] ! 

IGO = 4, Fourth-Order Runge-Kutta-Gill 

c = 1 - /2/5 c = -(2 + 3//2) 

1 o 

c 0 = -2 + 3/2 c c = 2 + /2 - 

2 6 

c 3 = 2 - /2 c y = 1/6 

c 4 = l + 1//2 c 8 = -1/3 
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f g i = f( yi’ t i ) 

| h 

L y = y L + 2 


' • § 2 = f( y» t i+i ) 

2 

2 s y = y + c 1 (§ 2 " s i )h 


g l = C 2 g l + C 3 g 2 


3 < 


g 2 = f (y.t i+ i) 

2 

y = y + c 4 (g 2 - g x )h 

§1 = C 5 g l + C 6 g 2 ■ 



g 2 = f( y’ t i+ i ) 


L y it i = y * <=782 + c 8 g i )h 

Figure 9 is a flow diagram of the numerical integration logic flow. 

Subroutine UPDATE 

Subroutine UPDATE computes the Earth ephemeris data and tracking station 
motion, provides the logic controlling computation of the observations and ob- 
servation partial derivatives, and computes the filtered state estimate. Figure 
10 shows a general logic flow for this routine. 

Subroutine MOTION 

In addition to calling subroutine ACCEL, which computes the dynamic accelera 
tions , the model compensation differential equations are evaluated. Further, the 
error covariance matrix differential equation is evaluated. Figure 11 is a flow 
diagram of the corresponding logic. See also Table 1. 


Subrout i no ACCEL 

. The true (simulated) and nominal accelerations (gravitational and thrusting) 
are computed according as the parameter, JP, equals 1 or 2 , respectively. This 
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includes both the gravitational and the thrust accelerations. For the nominal 
acceleration, the coordinate system used depends upon the value of the parameter, 

KSW . 

KSW = 1 , Heliocentric rectangular coordinates 

2 , Local rectangular coordinates 

3 , Local angular coordinates 
Figure 12 is a flow chart of the subroutine. 

Subroutine AMATRIX 

The partial derivatives of the state differential equations are computed. 

The forms of the derivatives depend on the coordinate system (KSW = 1, 2, or 3), 
and on the type of equations used to model the nominal thrust acceleration errors 
(MODOP - 11, 12, 21, 22). Figure 13 is a flow diagram for the subroutine.’ 
Subroutine HMATRIX 

The partial derivatives of the observations are computed and provided to 
subroutine UPDATE for use in the filter equations. Figure 14 is the flow diagram. 
Subroutine OBSERV 

The observations (range, range-rate, sun-planet angle, star-planet angle, 
and/or sun-star angle) are computed. Both the true and nominal observation values 
are determined. Figure 15 is the flow diagram. 

Subroutine GAUSS 

This routine provides random samples from a standardized normal (Gaussian), 
distribution. 

. Subroutine MINIPLT ' 

(See Department of Aerospace Engineering library documentation.) 




sin cos E 3 


T* + RCa'-’-'+c^] cos e 2 

sin e 2 sin Eg 

T* + sin ( e 2 + n 2^ COS ^ e 3 + ri 3^ 

R[a*+e 1 +n 1 ] cos (e 2 +n 2 ) 

sin (e 2 +n 2 ) sin (e 3 +n e ) 

sin cos Eg 
cos e 2 

sin Eg sin Eg 

sin ( e 2 +r| 2 ^ COS < ' e 3 +n 3^ 
cos (£ 2 +n 2 ) 

sin (E 2 +n 2 ) sin (Sg+Hg) 
for Vector Function, S 


T* + 

R[a*+e 1 +n 1 ] 


T* + R[a*+e 1 ] 




MODOP 


MODEL 


with IQ = 0 


Noise Mapping Matrix 


0 

0 ■ 

p = ap + pa t + rQr T 


E 

R 

I 



• rp 

r = 

1 

• 0 • 

• P = AP + PA + Q 



11 

1 



- 



k - T T 


COS ijj 

12 

4 

p = ap + pa + rQr 






R = 

sin ^ 

21 

2 or 3 








0 

22 

5 





Table lb. State Error, and State Noise Covariance 


I 

I 

-sin 
cos ij; 
0 


Logic 
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Code 

Dim 

Type 

| Text 

Description 

WEARTH 

' 1 

R 

W E 

Orbital angular velocity of the Earth 

AU 

i 

1 1 

R 

C 

1 Astronomical Unit (a.u.) 

PROBE 

18 

R 

X 

Initial value of state vector to be estimated 





(1) - (6): Position and Velocity, X Y Z X Y Z 

PROBER 

6 

i 

• R 


(7) - (24): e 1 e 2 e 3 n 2 n 3 a 2 a 3 ^ B 3 

Initial state vector errors in position and velocity 

GSUN 

| 1 

R 

M 

Gravitational parameter of the sun 

T ST ART 

1 ^ 

R 


Initial time 

TSTOP 

1 

R 


Final (problem termination) time 

TRANGE 

1 i 

R 


Range observation interval 

TRNGRT 

1 

R 


Range rate observation interval 

TSNPLT 

1 

R 


Sun-planet angle observation interval 

TSTRPLT 

1 

R ' 


Star-planet angle observation interval 

TSNSTR 

1 

R 


Sun-star angle observation interval 

THRSTAX 

1 

R 

a* 

Nominal value of thrust acceleration 

TI1KSTER 

i 

4 

R 


(1): Frequency (cycles/sec) of sinusoidal acceleration 

error variation 





(2): Percent of nominal thrust acceleration magnitude 

to be used as magnitude error 





(3): Standard deviation for random component of 

thrust acceleration error magnitude 



i 

1 

i 

j 

i (4): Standard deviation (radians) for random component 

of thrust acceleration pointing . error 

ACCELER 

' 6 ’ 

j 

R 

1 

i 

Not used 

PLATFER 

! 

i 6 1 

! 

R 

; 


Not used 

STERCOV 

i 

18 ! 
i 

i 

R 

i 

P 

Initial values of diagonal elements of error covariance 

STMSCOV 

! 

i 

, 18 1 
j 

i 

! 

R 

i 

] 

Q 

Constant values of diagonal elements of state noise 
covariance 

MOD OP ! 

— 1 

1 1 

i i 

] 

i j 
i 


Model option selection parameter 


Table 2 a. NAM Lb I DT/ 1 N PUTS / 
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Code 

! 

Dim 

Type 

Text 

Description 




i- 

i 

MODOP =0,1- Model 0 MODOP = 21 - Model 3 





MODOP =11 - Model 1 MODOP = 12 - Model 4. 

MSIMOF 

1 

I 


Not used 

OBSCOV 

5 

1 

R 


Nominal observation error covariance (elements 1-5 
correspond to the elements for TRANGE, TRNGRT, 

| TSNPLT, TSTRPLT , TSNSTR , respectively) 

TOBSCOV 

5 

R 


| True observation error covariance (elements cor- 
respond to those above (OBSCOV)) 

STAR 

3 

R 


Unit vector pointing to navigation star 

DTI NT. 

1 

R 

h 

Integration step size 

I PR NT 

5 

R 


Print frequency for each observation type 

UNMOD 

12 

R 


Initial values of model compensation parameters 

IPLOTYP 

i 

I 


Plot type flag: Values may be any number of digits 

from 1 to 7, in any order. Thus, 127 or 271, etc. 
results in plots on printer paper, teletype, or CRT, 
and microfilm. 

i 



! 

i 

1 

Codes are: 1 = line printer; 2 .= teletype; 3 = 

12" paper, ball point pen; .4 = 12" paper ink pen; 
5 = 30" paper ball point pen; 6 = 30" paper ink 
pen; 7 = 35 mm film- 

TOBS 


R 


initial times for each observation type 

TUP 

i 

R 

i 

. i 

Flag for controlling initial observation processing' 

TUP < 0. - Process first observation^ ) at initial 
time prior to integrating first step 
TUP > 0. - Integrate first step then process observations 

MAX IT 

• i 

• T 

i 

Not used 

E 

(3 ' 

R 

! 

E(l') contains age-weighting parameter, 1 < S < » , 

| if desired. E(2)-E(6) not used. 

1 : 

KAX 

1 

■; 

j I 

i 

| . ■ 

j ! Flag controlling coordinate system for compensating 

! ; models : 

i 


. 

i 

! 

1 : 

i 

j 

! • | 
j 

! 1 KSW = 1 Heliocentric rectangular coordinate .(X Y Z) 

■ ! KSW = 2 Local rectangular coordinates , (x y z) 

j : KSW - 3 Local angular coordinates (<J> y 0 ) 

L 

[ ‘ 

1 


Table 2a (Cont'd.) 



Code 

Dim 

Type 

Text 

Description 

IQ 



I 

Not used 

KAGE 

1 

I 


Flag controlling age-weighting suboptimal filter 
KAGE < 0, use suboptimal age-weighting (not used) 

PB 

1 

R 


On-line control parameter 
PB < 0, employ on-line display 

JPLOT 

1 

i 


Not used 

AUX 

6 

R 


Auxiliary acceleration error simulation parameters 

(1) = s (rate of motion of tip of accel . vector) 

(2) - y (max value of y specifying boundary) 

(3) = 0 (initial value of 0; updated at each y conta< 

(4) ’= (f> (initial value of <j>, direction of accel. 

(5) = a Q (standard deviation of the noise in 0) 

0 

(6) = a (standard deviation of . the noise In s) 

s 

BETA 

5 

R 


Not used 

DZT 

18 

R ' 


Not used 

i 

zs 

■ i 

3 

R 

ZS 

Z-component (geographic rectangular coordinates) of 
tracking station 

I 

1 

D 

3 

R 

. D. 
i 

/x 2 .+y 2 . 
si si 

W 

4 

R 

V+oi 

w(l) = Rotational angular velocity of Earth; 
w(2)-w(4) = Initial angular displacements of tracking 
stations 


Table 2a. (Cont ! d. ) 


Code 


Dim 


Type 


Text 


Description 


Complete 18 x 18 a priori state error covariance 
matrix (324 elements). If this NAMELIST is used 
NAMELIST/ INPUTS/ must precede /APCOV/; the array 
STERCOV must be zero, otherwise /APCOV/ will no!, 
be read 


Table 2b. NAMELIST/APCOV/ 



Table 4. COMMON/PARAMS/ 
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